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ABSTRACT 


An  analysis  of  a seal  model  is  made  where  the  rotating  element  has 
both  fixed  tilt  and  two-lobe  waviness.  The  stator  is  assumed  to  be  gimbal 
mounted  and  to  have  inertial  mass.  Hydrodynamic  lubrication  is  assumed, 
following  the  short  bearing  or  narrow  seal  model.  Conditions  are  examined 
where  the  stator  precesses  in  synchronism  with  the  rotor  rotation.  Par- 
ticular interest  is  given  to  operating  conditions  where  such  behavior  appears 
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INTRODUCTION 


In  recent  years,  numerous  theories  have  been  advanced  to  investigate 
the  mechanism  of  mechanical  face  seal.  It  is  now  widely  accepted  that 
misalignment  alone  or  coupled  with  seal  surface  waviness  is  the  cause 
for  the  existence  of  a lubricating  fluid  film  between  seal  faces.  In 
the  case  of  non-flexible  mounted  face  seal,  while  stable  operation  is 
possible  if  only  one  of  the  surfaces  is  misaligned,  Harrdt  & Godet  [l] 
have  shown  that  axial  vibration  may  occur  if  both  of  the  faces  are  mis- 
aligned. Flexibly  mounted  face  seals  are  used  in  many  applications  to 
accomodate  misalignment.  These  should  have  at  least  second  mode 
wavlness  on  one  of  the  surfaces  in  order  to  have  stable  motion  has  been 
stated  and  experimentally  proved  by  S tanghan-Batch  & Iny  [2l. 

Thin  film  flows  between  solid  boundaries  can  cause  viscous  heating 
which  may  lead  to  thermal  deformation  of  the  solid  surface  and  consequent 
alterations  in  the  flow.  The  heating  and  deformation  may  enhance  each 
other  progressively  to  cause  large  surface  distortions  and  failure.  The 
phenomenon,  referred  to  as  "Thermoelastic  Instability"  has  been  studied 
by  Banerjee  [3"l.  Critical  sliding  speed,  above  which  instability  may 
occur,  has  been  found,  for  thermal  conductor  sliding  on  insulator: 


V = h <l-~- 
crit 


(1) 


where  mean  f i lm- th i ckness  h is  held  fixed,  K is  the  wave  number  of  the  dis- 
turbance, K the  thermal  conductivity  of  the  conductor  and  a its  coefficient 
of  thermal  expansion.  The  critical  speeds  measured  experimentally  were 
found  to  match  those  predicted  by  Eq . (1)  extremely  well.  This  equation  is, 
however,  restricted  to  constant  h and  non-flexible  mounting.  Later,  an 
improved  model  where  constant  h condition  is  relaxed  has  been  studied  [4], 
finding  that  thermoelastic  effects  lead  to  a continuous  change  in  surface 
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waviness  as  well  as  mean  f i lm- th ickness , with  changing  speed.  Conse- 
quently, the  unstable  region  cannot  be  obtained  in  practice  because  there 
exists  no  path  to  it  from  the  initial  waviness  at  rest.  This  theory  is 
based  on  the  assumption  that  mounting  is  inertialess  therefore  eliminates 
the  effect  of  misalignment  completely.  In  a set  of  experiments 
later  performed  [5],  where  a gimbal  has  been  used  to  simulate  the  flexible 
mounting,  although  mean  film-thickness  is  predict’  i very  well  by  the  theory, 
it  is  observed  that  first  mode  wave  grows  almost  three-fold  in  the  tested 
speed  range.  This  increase  is  believed  to  be  thermoelastic  in  nature, 
and  is  a consequence  of  gimbal  inertia,  loading  the  faces. 

In  view  of  the  above  fact,  a new  model  for  the  flexible  mounted  face 
seal  is  studied  here,  taking  both  first  and  second  mode  waviness  into 
account,  to  provide  a better  understanding  of  its  dynamic  characteristic 
and  thermoelastic  effect. 


METHOD  OK  APPROACH 


A three  degree  of  freedom  fact*  seal  model  has  been  proposed  which 
is  illustrated  in  Fig.  (1).  The  "rotor"  is  a floating  element  with  both 
tirst  and  second  mode  waviness  on  the  seal  surface;  the  "stator"  is  tlex- 
iblv  mounted  such  that  the  edge  of  its  flat  surface  can  have  angular  dis- 
placement about  its  center  of  mass.  It  is  also  assumed  that  rotor  is  a good 
thermal  conductor  while  stator  is  a good  insulator.  litis  assumption  is 
reasonable  [hi  and  leads  to  considerable  simplification  of  the  analysis. 

lhe  procedure  for  solving  the  problem  can  then  be  stated  as  follows: 
[refer  to  block  diagram.  Fig.  (2)1 

1.  Assuming  that  film  stiffness  largely  insures  tracking,  that  is,  the 
stator  will  follow  the  rotor  synchronously  (but  not  necessarily  in 
phase)  . 

2.  Considering  steady  state  solutions  only  and  choosing  the  most  con- 
venient rotating  coordinates. 

3.  Assuming  small-perturbations,  thus  allowing  hydrodynamic  pressure 
to  be  a linear  function  of  operating  waviness. 

Case  1 . Removing  second  mode  waviness  from  rotor. 

Tli  is  enables  one  to  calculate  the  force  and  moment  integrals. 

This  is  examined  to  see  what  is  the  role  of  the  first  mode  film- 
thickness  when  the  stator  is  flexibly  mounted. 

Case  2.  The  configuration  where  first  and  second  modes  exist. 

The  Information  previously  obtained  is  used  to  simplify  the 
analysis.  It  will  be  found  that  the  problem  can  still  be  treated 
as  a linear  one  where  force  and  moment  depend  "seperately"  on  second 
mode  and  first  mode  operating  waviness,  so  long  as  waviness  amplitude 
is  small  compared  with  film-thickness.  Once  force  and  moment  integrals 
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have  been  carried  out,  one  can  easily  arrive  at  the  last  step  of  the 
block  diagram,  where  operating  waves  are  composed  of  initial  waviness 
and  thermoelastic  deformation  on  the  rotor,  plus  stator  tilt  which  is 
caused  by  hydrodynamic  forces.  Use  of  this  result  enables  one  to  obtain 
the  condition  of  thermoelastic  instability,  and  to  find  out  critical 
sliding  speed. 

When  the  restriction  of  smal 1 -perturbat ion  is  removed,  it  will  be 
shown  that  the  characterist  ics  of  the  problem  remain  the  same  while  only 
a numerical  correction  factor  need  to  be  considered. 

Reiterating,  the  motion  of  face  seal  elements  will  be  investigated, 
where  the  rotor  is  rigidly  fixed  to  a shaft,  and  may  be  tilted  as  well  as 
wavy  in  the  second  mode.  It  will  be  treated  as  a thermal  conductor  and 
subject  to  thermal  deformation  as  the  result  of  frictional  heating  by  the 
fluid  film.  The  stator  will  be  treated  as  a flat,  thermal  insulator 
supported  in  gimbals.  Operation  will  be  studied  where  the  waviness  component 
of  f ilm-thickness  is  small  relative  to  the  mean  value  (a  condition  observed 
in  experiments) , and  particular  interest  will  be  given  to  those  operating 
conditions  where  this  configuration  becomes  impossible.  Hydrostatic  pressure 
and  leakage  effects  will  be  omitted  from  the  present  analysis. 

Comment  on  Notation.  In  that  which  follows  the  tilda  (~)  will  denote 
quantities  associated  with  the  stator  and  the  carat  (~)  will  denote 
quantities  associated  with  the  rotor.  First  mode  or  tilt  of  the  stator 
would  be,  in  such  notation,  x cos  9,  when  0 is  a measure  of  angular  position. 


IUKORY 


Hvdrodvnamic  equal  ions  : 


Referring  to  Fin.  (3),  it  will  bo  assumed  that  radius  K is  much 
larger  than  width  L ot  seal  surface,  therefore  a one  dimensional  Reynold's 
equation  with  narrow-face  approximation  [7]  can  he  written  as 

•£-  (h  3 ^1  - b + 12  (21 

where  h is  overall  fi lm- thickness , a function  of  x and  t but  not  of  v. 
Hence,  integrating  p twice  with  boundary  conditions  p « 0 at  v - fL/2 


results  in 


p - C-Lv2  . J-l2)  ^ an  + & 

1 ' 2 ' B \]  9x  ,3  9t 


Width-averaged  pressure  can  be  calculated  as 


P - - pdv 


3 dx  ’ . 3 St 
2h  h 


Rotating  coordinates  mav  be  chosen  to  simplify  the  problem.  Under 
the  assumption  of  synchronous  motion,  it  can  he  recognized  that  the  pressure 
distribution  is  rotating  wave,  but  non-varying  in  shape;  therefore  it  will 
be  convenient  to  choose  moving  coordinates  such  that  pressure  function  will 
be  fixed  with  respect  to  them.  In  this  favorable  reference  frame  system, 
the  film-thickness  wave  is  fixed.  This  does  not  require  that  the  wave 
representing  tilt  of  the  rotor  or  stator  to  bo  in  phase  with  the  film- 
thickness,  but  only  that  the  phase  relationships  in  the  ensemble  are 
preserved  in  time.  Recall  that  the  sur face  waves  are  fixed  to  the  rotor, 
hence  dh/St  disappears,  and  the  stator  surface  is  moving  past  with  the 
velocity  -V.  When  this  sign  is  taken  into  account  Kq.  (41  is  replaced  by 

pVl,2  *h 

P - o-  (51 

2h3 
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where  V is  now  simply  the  magnitude  of  the  sliding  speed.  This  equation 
shows  that  pressure  is  a function  of  position  and  that  squeeze  effect  is 
not  seen  in  the  chosen  coordinates. 

Force  and  moment,  according  to  definitions,  are 

r - 

f = J+  pLRde  (6) 

P 

My  = -J+  pLR2sin6d0  (7) 

P 

My  = J+  pLR^cosQdG  (8) 

P 

where  "p+"  is  used  to  denote  that  integrals  are  effective  only  when  pressure 
is  positive,  while  in  the  region  of  negative  pressure  cavitation  is  assumed 
to  occur  and  pressure  is  zero.  This  condition  has  been  experimentally 
pioved  to  be  quite  accurate  for  the  face  seal  [2], 


Motion  ot  Ht.ilor  :uul  Jo  1 1 n 1 1 1 on  o l_  lot  a t inn  coord  I nat  oh 


Ulion  seal  Is  .steadily  operating  at  certain  speed,  hydrodynamic 
torro  Is  const. snt  In  amplitude  but  moving  along  tbo  edge  ol  the  stator 
surface.  It  can  bo  seen  that  it  the  tilting  amplitude  t.  (Kin.  4)  is 
small  compared  with  K,  which  is  Rent' rally  true,  then  the  angle  ol  tilt 
will  he  approximately 


Y 


('» 


Rotation  coordinates  ( see  Kin.  •«)  mav  l>o  chosen  so  that  the  v axis 
is  the  axis  about  which  tilt  occurs,  and  it  mav  also  he  taken  to  lie  In 
the  plane  ot  the  lace  ot  the  stator.  the  x axis  is  porpendlculat  to  the 
v axis  and  points  in  the  direction  ol  maximum  elevation  ot  the  stator.  flic 
axis  is  that  about  which  the  rotor  turns  and  is  taken  t be  vertical, 
while  the  t.  axis  represents  the  normal  to  the  plane  ot  the  stator  lace. 

llie  dynamics  ot  an  axially  symmetric  mass  movinn  with  such  an  K.ulcrtan 
triune  are  well  understood  (see  Appendix  K lor  further  comment).  Kor  the 
case  where  the  tilt  annle  Y (measured  between  the  e and  /.  axes)  retains 
constant  magnitude,  while  the  direction  ol  the  v axis  rotates  steadily  at 
tin'  angular  speed  ot  the  shaft  , w,  the  moments  acting  on  the  mass  are  such 
t hat 


M - 0.  M - 0 
x * 

M “ I u'*slnycosY  (III) 


It  the  stator  is  taken  to  be  a thin  ring 

l - mK‘/2  ill) 

where  K Is  the  moan  radius  ol  the  ring,  and  Kq . ill')  become# , lor  small  y 


M 

v 


m V* 

~l  R 


(l*) 


Recall  that  this  moment  is  ot  constant  magnitude  and  always  acts  t o oppose 
the  angular  displacement  y. 
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Heating  ettect  a ml  the  nnoe  last  tc  detormatlon 

1(  has  been  shown  ["  8 ] that  for  I ncompre s s l b le  1 low  in  thin  films, 
viscous  heat  generated  within  the  lilm  Is  given  bv 

q - J dz 

o 

Under  tin'  short  hearing  model,  the  fluid  velocity  distribution  in  the 
film  approaches  Couette  flow  which  is 


u 


v£ 


therefore 


q 


(1  ') 


Georgopoulos  has  found  that  under  typical  seal  conditions  converted 
heat  is  only  a small  amount  of  percentage  of  that  conducted.  In  view  ol 
this  It  can  he  assumed  that  generated  heat  will  he  removed  bv  conduction 
alone,  furthermore,  it  would  conduct  totally  Into  the  rotor  which  is  the 
only  conductor  under  consideration. 

The  equation  relating  heat  input  to  a surface  and  the  corresponding 
surface  curvature  has  been  shown  bv  Hurton  et  al  I 10 1 , and  will  be  written 


as 


> 

dx* 


K 


(I/O 


where  $ , means  the  thermoelast tc  deformation  on  the  rotor  surface,  and 
th 

q*  is  non-uni  form  part  of  total  heating  q since  uniform  heating  can  only 
lead  to  zero  curvature'. 
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ANALYSIS  KOK  SMALL  - PKRTlfRHATI  (>N 
Case  1 . Rotor  has  only  tirst  mode  wavtnoss 

In  tills  case,  wi'  shall  examine  who t her  a s toady  stato  solution  Is 
posstblo  wlion  thoro  exists  only  I list  modo  wavtnoss,  and  Its  amplitude 
vory  small  compared  with  moan  t'i  lm- thickness  li . 


Overall  t l lm- th Icknoss  Is  the  sum  of  h and  the  operating  wave  li  ^ , 


where  h comes  I rom  rotor  wavtnoss  and  stator  motion.  To  he  consistent 
o 


with  chosen  coordinates,  it  will  he  written  as 


ti  - e sln(&H>  ) - z cosO 

o 1 11 


(IS) 


where  z j denotes  amplitude  ol  rotor  wavtnoss,  / ^ amplitude  ot  the  stator 


and  Is  a phase  angle  determined  hv  initial  waviness  and  t liermoe last  1 c 


deformation.  Kor  convenience,  h ^ may  he  simplttled  Into  one  term  as 


h « z ' s 1 n t (H  ♦ , 1 
o 1 1 


(lh) 


with  the  relation 


zj  - (zj  - 2*^  sin  ‘ 


(17) 


and 


Z.jslll'lj  “ ZjSiu^j  - Zj 


(18) 


Recalling  Kq . (S)  and  noting  that  positive  pressure  implies  dh/dfl  > 0,  or 

(19) 


at  j cos  (&♦■•)  > 0 


'therefore  the  positive  pressure  region  around  the  seal  surface  will  he 


- 2 - ♦ 

2 l 


- - * 

2 1 


With  use  ot  this  Interval  and  the  assumption  ot  small  amplitude  waviness 
(such  that  *j/h  « 1),  the  following  results  can  he  easily  derived  I rom 


Kq.  (h) , (7)  and  (8) 
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z-smr 


—3  1 ' 

4hJ  1 


n^VL  R , 

\ = -T3-  Zjcos^  (22) 

4h 

Remembering  = 0 and  My  > 0 for  the  chosen  coordinates,  it  follows  that 
must  be  zero.  With  this  constraint,  Eq.  (18)  becomes 

Zj  - z^sinijj  (23) 

inserting  this  into  Eq.  (17)  yields 


Z1  ” V* 1 * *1  sin6 


Physically  = 0 means  the  stator  will  follow  (tilt)  in  a way  such  that 
the  operating  wave  h^  is  90°  out  of  phase  with  the  stator  wave.  The 
generated  pressure  will  be  symmetric  about  X axis,  therefore  it  will 
produce  no  moment  about  this  axis.  Also  it  is  observed  from  Eq.  (23) 
that  z^  should  be  smaller  than  z^.  Consider  the  special  situation,  where 
the  stator  is  fixed,  then  the  operating  wave  reduces  to 

h = z sin9  (25) 

o 1 

which  is  the  case  investigated  earlier  [4].  Another  limiting  case  happens 
when  Zj  ® z^,  it  turns  out  that  0^  = tt/2  and  zj  = 0 which  means  stator 
keeps  parallel  with  rotor,  obviously  no  pressure  will  be  generated  in  this 
position  hence  impossible.  Nevertheless  it  points  out  that  if  zj  is  very 
small  compared  with  z^,  then  z^  approaches  and  there  is  a slight  phase 
shift  from  0^  * tt/2. 

The  relation  between  force  and  moment  derived  from  Eq.  (20)  and 


Eq.  (22)  is 
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inserting  Into  Kq.  (12)  yields 


nR 
2 mV  2 


(2b) 


Rewrite  this  by  changing  m to  w/g,  and  ^ ■ N,  where  w is  the  weight  ot  the 


stator,  N is  revolutions  per  minute,  and  it  becomes,  after  some  rearrangement 


‘i  - <£>  6 <T>2  (27) 

w 

In  general  the  axial  load  f is  not  necessarily  related  to  stator  weight. 

In  this  case,  however,  we  may  assume  that  axial  load  is  equal  to  the  weight 
of  the  rotor  and  be  consistent  with  the  experiment  previously  performed. 

It  is  then  found  that  z^  will  be  too  large  even  for  relatively  high  speed. 

For  example,  take  f = lOw,  N *=  3000  r.p.m.,  tilt  amplitude  will  be  approx- 
imately 

z i = 1 . 5mm 

and  for  N » 1000  r.p.m., 

Zj  “ 13.75mm 

But  as  stated  before,  steady  operation  requires  z^  always  bigger  than  'z 
which  means  initial  waviness  should  have  at  least  same  order  of  amplitude 
as  z.j,  and  is  impractical. 

It  can  be  concluded  that  for  a flexible  mounted  face  seal,  if  axial 
force  is  comparable  with  stator  weight,  then  there  will  be  no  steady  state 
solution  under  the  assumption  ot  synchronous  motion.  The  reason  lies  in 
the  fact  that  single  wave  provide  both  force  and  moment.  Hie  moment 
required  by  the  stator  is  very  small  while  the  moment  generated  is  large  ii  same 
wave  has  to  support  the  load. 

Returning  to  F.q . (27)  and  requiring  'z^  * 6^ , one  finds  the  limiting 
speed  below  which  operation  is  strictly  forbidden,  this  being 

n2  . .Mi2  _[fi_ 

w 8tt  — 
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Case  2.  Rotor  has  both  first  and  second  mode  waviness 

The  operating  film  wave  in  this  case  can  be  expressed  as 

hQ  - ZjSin(&f0j)  + z2sin(2&f®2)  ' ZjCosG  (29) 

where  z 2 is  the  amplitude  of  second  mode  on  rotor,  ®2  is  the  phase  angle 

relative  to  the  stator.  As  in  previous  cases,  it  can  be  reduced  to 

h = z'sin(frf$  ) + z sin(26+0  ) (30) 

o 1 12  2 

where  zj  and  are  defined  in  Eq.  (17)  and  (18).  The  condition  of  positive 
pressure  becomes 

^ = z jcos (8+$j)  + 2z2cos(2&4-tf>2)  > 0 (31) 

— 3—3 

Assuming  again  h^/h  <<:  1,  amd  1/h  1/h  , force  and  moment  equations 

assume  the  following  forms 

3 ■ 

f «*  j+  [zjcos(frf$^)  + 2z2cos(2&f02)  ld9  (32) 

2h  p 

3 

^ **  ^^3  R J+[Z|Cos(&f<1)  + 2z2cos(2&f(#2)  IsinGdG  (33) 

2h  p 

3 

My  =*  ^^3**  |+  [zjCOs((H-$j)  + 2z2cos (26+®2)  IcosGdG  (34) 

2h  p 

Eq.  (31)  shows  that  integration  interval  p+  is  a function  of  z j , z2> 
and  ®2 , therefore  the  analytical  treatment  of  Eq . (32),  (33)  and  (34)  seems 
impossible.  However,  conclusions  from  the  previous  case  and  the  symmetric 
property  of  second  mode  can  be  used  to  make  reasonable  simplifying  assumption 
which  will  make  closed  form  solution  possible. 

As  stated  before,  if  first  mode  alone  is  responsible  for  the  face  load, 
then  the  moment  generated  will  usually  be  too  much  for  equilibrium.  In 
order  that  moment  be  of  reasonable  magnitude,  one  must  let  zj  become  very 
small  such  that  it  is  unable  to  lift  the  load.  On  the  other  hand,  the  second 
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mode  has  the  ability  to  lift  load  while  producing  no  moment  owing  to  its 
symmetric  property.  So  in  the  combined  case  of  first  mode  plus  second 
mode,  it  is  reasonable  to  assume  that  zj  should  be  very  much  smaller  than 
second  mode  waviness  , such  that  the  second  mode  is  largely  responsible 
for  the  load  support  and  the  first  mode  wave  is  small  but  sufficient  for 
dynamic  tracking. 

To  prove  the  above  argument,  recalling  the  dynamic  equation  of  the 
stator  again,  it  becomes,  for  this  case 


~ r2  l3  » z' 

Zj  = z2  J+[-j-cos(&f®1)  + 2cos(2{H-$2)  ]cos0d9  (35) 

mV  h p 2 

where  and  0^  must  satisfy  condition  = 0.  By  defining 

zi 

C = r1  (36) 

z2 


and  substituting  Eq.  (32)  into  Eq.  (35)  it  will  result  in,  after  some 
rearrangement, 


f o <[+  [£cos(6+$  ) + 2cos(20+0„) 1cos0d8 

z = (-“-)  (-)  (— ) £ - ±— 

1 "-2  ~ N j+  [Ccos(frf$1)  + 2cos(2&f02)]d9 


2 TT  w 


P 

. (60)2  INTfla. 

^o_2''v~mN;  4nt(f); 


(37) 


2tt  w 

where  INT(M)  and  INT(F)  are  notations  for  moment  and  force  integrals. 
Computer  results  [Fig.  (5)  & Appendix  Al  indeed  show  that  Z must  be  much 
smaller  than  unity  to  keep  z^  within  practical  range.  For  instance,  if 
N = 3000  r.p.m.  and  f/w  = 10,  then 


*.2 

'‘INT(F)  ' 


(cm) 


for  Zj  S 2.5  x 10-'  fc»)  ,INT(M) /INT(F)  should  be  less  than  1.25*10’2.  From 
Fig.  (5)  one  finds  that  C < 0.1  (for  any  possible  phase  angle)  or  zj  < 0.1 
And  for  N ■ 1000  r.p.m.,  Z should  be  less  than  0.01. 


N> 
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As  shown  before,  small  £ which  implies  small  zj  (»  Cz2)  will  mean 
that  the  stator  follows  in  the  position  almost  parallel  to  the  rotor. 

It  can  be  expected  that  if  the  stator  approaches  the  inertialess  condition 
then  zj  approaches  zero  and  second  mode  will  be  the  only  waviness  in  the 
operating  film. 

Under  the  assumption  of  small  Q,  positive  region  can  be  determined 
approximately  by  only  second  mode  since  first  mode  can  only  slightly  change  it. 
Hence  p+  will  be 

TT  ^2  - TT  *2  3tt  02  5TT  ®2 

*4'T  9 4’T  dT'T  e T‘T 

Assuming  that  h^/h  « 1,  and  £ « 1 , force  will  be  approximately 


3 

f = z2  J+  2cos(29+02)d6 


2h  p 

2u,VL3  « 


h3 


(38) 


this  equation  fits  very  well  with  numerical  data  [Fig.  (6)  and  Appendix  A 
with  e = OJ.  Mean  film-thickness  can  be  derived  from  it  as 


h - (aaL.) 1/3  (s2) 1 


/ 3 


(39) 


^ !+  z|cos(Uf$^)  sin9d9  + [*+  2z2cos(20f02) sinOdoj  (40) 


With  same  kind  of  reasoning,  the  moment  can  also  be  calculated.  From  Eq.  (33) 
3„ 

2h  p p 

The  second  integral  can  be  equated  directly  to  zero  by  recognizing  the 
symmetrical  property  of  the  second  mode,  the  result  is 

Mx='^-3"  r 
x *■  «•  *■ 

For  the  chosen  coordinates,  * 0,  hence  there  exists  a relationship  between 
$1  and  tf2 , which  can  be  derived  from  Eq . (41)  as 


(41) 
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with  the  above  restriction,  M^,  will  he  the  total  moment,  which  it  follows 


the  same  procedures  as  did 


» ^VL-  R z'  'cos ($-£„)  + tcos*  ] (43) 

^ 2h3  1 L 1 2 2 II 

I 2 ^ 

Computer  results  [Fig.  (7)  and  Appendix  Al  show  that  the  total  moment  is 

accurate  even  when  £ approaches  0.5.  From  Kq.  (12),  tilt  of  the  stator 


becomes 


2mVh 3 1 ' 


z'  ' cos($ ,-0»)  + -z  cos$  f 


Equation  (39)  and  (44)  show  that  mean  f i lm- thickness  depends  on  second 
mode  only,  and  tilt  of  the  stator  depends  on  first  mode  wave  amplitude 
alone  but  will  he  affected  by  second  mode  position  (tf^)  • 

Heating  effects  and  thermoelastic  deformations  will  now  be  dealt 
with.  Recalling  that 


when  h h , it  will  become  approximately 

, - *o  (45) 

h(l-H^)  h h 
h 

where  first  term  is  the  uniform  heating  which  does  not  contribute  to  surface 
curvature.  The  second  term  is  non-uniform  heating  caused  by  operating  film 
wave.  It  can  be  further  expressed  as 


q'  ■ |zjsin(&f*j)  + z2sin(2&t<2>2) 
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Subs  t i tut  ing  into  Eq . (14)  yields 


d2$ 


2 

r1  = |z|sin(6+$j)  + z,sin(2&4-02>) 


d x“  Kh 

After  integrating  twice, 

2„2 


(47) 


Sth  - |z'sin(&4-$1)  + -^sin(2&f«2)|  + CjG  + C2 

Kh 


(48) 


Clearly,  both  and  C2  should  be  zero  because  of  seal  geometry. 
Therefore  thermoelastic  deformation  is  composed  of  first  mode  and  second 
mode  component  which  will  be  expressed  as 


.2  2 


*th.  = zisin(&fV 

1 Kh 


’ , = -4  z sin(2&f0  ) 

th2  4Kh2  2 2 


(49) 


(50) 


It  is  noted  that  the  above  derivation  implicitly  assumes  that  heating  is 
generated  by  the  full  film.  Stanghan-Batch  [2]  has  shown  that  although 
cavitated  regions  exist,  the  fluid  is  partially  continuous  around  seal 
surface.  Therefore  actual  values  of  heating  and  surface  deformation  should 
lie  between  the  extremes,  full  film  and  absent  film  when  p = 0. 

Referring  to  the  block  diagram  again  [Fig.  (2)1,  two  equations  can 
be  written  from  initial  and  operating  conditions  by  separating  first  and 
second  mode  components 

(51) 


z.  sin(9+^f  ) -f  S - z.cosQ  = z'sin(£H-$,) 
it  1 thjl  X 1 

z2isin(2frH2)  + Sth  = S2sin(2&4^2) 


(52) 


where  z^,  are  initial  amplitude  of  first,  second  mode  waviness  of  the 
rotor,  and  t|r 2 are  phase  angles  relative  to  the  stator  wave.  Substitute 
Eq.  (43,  49,  50)  into  Eq . (51)  and  Eq . (52)  which  yield 


z1}sin(&4-ijr  - zj  ^1  - ^ — j8in(£M-$j)  + z^cos© 

Kh 


l mVh3 


f,  uaV2R2  ' 


zi\  2 1 L " — 

F(*1,®9)viR  L Kh 


2~3  L1  ~=2“jsin(e+V  + COS0i 


whore 


2 2 

i in  (29-Hjr .,)  * z9  ^ s in  (26+09) 


F(#i 3 cos($1'®2)  + ^ cos$1 


From  Kq.  (54)  it  is  found  that 


+ 2 = 02  F 2tt 


which  moans  thermoolastic  deformation  in  the  second  mode  will  be  in  phase 
with  the  initial  waviness.  Operating  waviness  z2  can  be  obtained  by 
substituting  Eq . (39)  into  Eq.  (54)  as 

22  3-2/3 

uaV  R ,2uVL\  ^ ,1/3 

Z2i  ' Z2  * 4K  { f ) (z2>  (56) 

Therefore  if  initial  amplitude  z.,^  is  known,  then  operating  amplitude  z9 
can  be  determined  and  will  depend  on  load,  sliding  speed  and  other  physical 
propert ies . 

The  initial  conditions  of  the  rotor  waviness  including  amplitude  and 
phase  relationship.  Phase  shift  between  first  mode  and  second  mode 
wav i ness  will  be  defined  bv 


2 2 

^i  ’ T ' *1  2 ‘ *1  * n 


where  Atjr  M ranges  from  -tt/2  to  n/2. 
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Inserting  this  into  Eq . (53)  and  taking  two  orthogonal  components  yields 
two  equations 


0 


'li 


2 ~ r 

COsC-ptTT-A^)  = Zjj 


—3 

2mVh 


HR2L3F($1,02) 


(1-  ) ;COS$ 

Kh 


(58) 


0. 


- , /2,_  x ~ ~ T 2mVh3 

Zu  sln(TfcTT'^i)  ‘ zi  = zil'XI 


(1  - MW  R ) ,sin$  (59) 
M.R“LJF($,  ,0.)  Kh 


1 '2 

where  02  ranges  from  -tt  to  TT.  Recalling  Eq.  (42)  which  relates  to  02> 
and  dividing  Eq.  (59)  by  Eq.  (58)  will  yield 

0, 


~ tt  ^2  *2 

zi  jsin(-y±TT-A \|r.)  + s i n (-y±rr+A^  ) 


'li 


^ - COS02 


By  defining 


A = 


2mVh3 

_ 2 3 
URL 


Kh  crit 

and  substituting  Eq.  (60)  into  Eq.  (58),  one  finds 


02  _ 

cos(-y±TT-Ati)  (j-COS02)F(§]  ,02) 


...  0~ 

[-y  8in("y±TT-At  ^)  + Sinf-yiTT+A^)  ]cos$j 
Also  Eq.  (44)  now  becomes 

F 


= AB 


Z1  A Z1 


therefore 


AZU  ZliF<W 


(60) 


(61) 


(62) 


(63) 


(64) 


(65) 
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A computer  program  [Appendix  Dl  has  been  written  to  evaluate  AB  and  zj/Az^. 
Part  of  the  results,  where  negative  z^/z^  is  not  allowed,  are  shown  in  Table 
1.  It  is  found  that  AB  may  have  positive  and  negative  values  for  any  initial 

/v  a O 

phase  angles,  and  is  zero  when  * z^  and  ^ * 90  , which  means  zj  is 
restricted.  Although  this  is  true,  it  is  also  found  that  zj/Az^.  reaches 
maximum  (when  A^  = 0°)  or  near  maximum  (for  other  At  . ) at  AB  = 0.  Values 
of  z|/Az^  for  AB  = 0 are  shown  in  Fig.  (I).  Whether  these  values  can  be 
reached  depends  on  A.  For  example,  take  the  worst  case  when  At ^ = 0°,  Fig. 


(8)  shows  that  zj/Az^  = 1.75,  so 


~ -3 

z ' mV  . h 

— = (1.75A)  = 3.5 


'li 


B=0 


2 3 
p,R  L 


(66) 


By  substituting  for  h from  Eq . (39), 

z\  7mV2 


'li 


fR 


crit  * 

o z 

2 2 


(6  7) 


A limiting  case  is  that  m becomes  very  large  (approaching  the  fixed  mounting 
case),  then  zj^  will  approach  infinity,  which  is  the  case  investigated 
earlier  [l].  Eq.  (67)  can  be  changed  to 

.2  z! 


crit 


fR 


~ A A 

7m  znz2 


(68) 


Since  the  stable  condition  requires  zj  ^ z it  is  instructive  to  choose 
values  of  zj/z2  as  the  limiting  condition  to  evaluate  critical  speed.  For 
instance,  if  z{^z2  = 0.1,  then  Eq.  (68)  becomes 

.2 


fR 


crit 


70mz 


li 


The  ranges  of  f,R,m  and  z^^  may  vary  broadly  for  face  seals,  and  so  will 
the  critical  sliding  speed.  Eq . (69)  points  out  that  bigger  f,R  and  smaller 
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m.z^  will  make  ^cv^t  higher,  hence  safer,  and  vice  versa.  Also,  by 
noting  that  larger  f/m  requires  smaller  z £ relative  to  it  follows  from 
Eq.  (68)  that  the  radius  of  the  stator  and  the  initial  amplitude  of  the 
rotor  are  the  most  important  factors  in  considering  thermoelastic  effect. 

It  is  noted  again  that  when  m approaches  zero,  there  exists  no  operating  first 
mode  wave  in  the  film;  same  result  will  occur  when  initial  first  mode 
waviness  approaches  zero.  In  both  cases,  there  will  be  no  thermoelastic 
deformation  in  the  rotor  first  mode  hence  no  thermoelastic  instability. 

These  observations  can  be  reflected  from  Eq.  (69)  by  noting  that  Vcrit 
approaches  infinity  and  are  consistent  with  the  previous  theory  [4], 


RKSULTS  AND  DISCUSSION 


For  flexibly  mounted  face  seals,  if  there  exists  only  first  mode 
waviness  on  the  rotor,  and  axial  load  is  comparable  with  stator  weight, 
then  there  will  be  no  admissible  steady  operation  under  the  assumption 
of  synchronous  motion.  If  there  exists  only  second  modi-  waviness,  then 
thermoelastic  instability  will  not  occur  [4].  For  the  case  where  both 
first  and  second  mode  wavlness  Is  present,  the  condition  for  steady 
operation  requires  zj  ^ • It  is  found  that  axial  load  is  supported 

largely  by  the  second  mode  while  the  first  mode  wave  will  be  responsible 
Cor  the  motion  of  the  stator.  In  such  situations,  thermoelastic  deformat ion 
will  grow  when  sliding  speed  increases.  Although  It  *111  not  roach  infinity 
as  in  the  case  of  the  fixed  mounted  face  seal  [l],  it  may  reach  a certain 
value  where  the  requirement  of  z'  « fails.  This  enables  one  to  determine 
the  critical  sliding  speed,  which  will  depend  on  load,  mass  of  stator,  radius 
of  stator,  initial  first  mode  wavlness  of  the  rotor  and  initial  phase  angle 
between  first  and  second  mode  wavlness  of  the  rotor. 
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Figure  2.  Block  Diagram  for  Tliermoelastic  Analysis 


L 


INT(M) 

INT(F) 


X «■  0 


X = tt/6 


X - tt/3 


X = n/2 


amplitude  of  first  mode  wave  (zp 
amplitude  of  second  mode  wave  (Z2) 


Figure  5.  (Moment  integral/Force  Integral)  vs.  C for 
small-perturbation  Analysis 


1 


Figure  8.  ii/As^  va.  i at  AB  ■ 0 


Appendix  A - Program  for  computing  force  and  moment  Integrals 


PROGRAM  MOMEN  ( I NPUT » OUT  PUT ) 

R£  AL-fU,  > rtY  1 » H Y 2 » t AKOA 

C EPS1  - 2ND  MODE  / MEAN  FILM  THICKNESS 

C ZETA  - 1ST  MODE  / 2ND  MODE 

PI-3. 1*15926 

RE  AO — l.D.ZETA.EPSI 

1 FORMAT ( F6.0#F10.8«F3. 1 ) 

OTH-PI/D 

00-  LOO — I<U10 - 

LAMOA-O. 

K-0 

PRINT  2,  D.ZETA.EPSI 

2 . FORMAT ( />1X^F6.0»F ia.8.F 10.2 1 - - 

3 THETA--OTH 
MX  1 -0 • 

MX2-0. 

MT1-0. 

MY2-0. 

Fl-O. 

F2-0.  — - - 


THETA-THETA+OTH 

P-ZETA*COS( THE T A* L AMD A » ♦ 2 • OS ( 2.* THETA ) 

IF  ( P.LU-4. .AKO-. — THEJ-A-  .4.U  2.*PIJ  - GO -TO 


50 


100 


OF  1- P *0TH 
F1-F1+0F1 

0MX1«P*SIN(  THETA-l*OTH  - - - - 

0MY1-P*C0S( THcTA)*DTH 
MX  1 -MX 1+  OMX 1 
MY  1-MY1+DHY 1 

H-l.Vz  EtX*S1  nTThE  TA  ♦ L A MO  aT+E  PSI*  S I N(  2 . ♦THETA  ) 
HCUBE -H*H*H 
OF  2-OF 1/HCUBE 

F2»F2*DFZ-  - - - - 

0MX2-  DMX1/  HCUBE 
0MY2-  0MY1/  HCUBE 
MX2-MX2+0MX2 

-MY2-MY2+-0M-W 

IF  ( THETA  .IT.  2.*PI)  GO  TO  5 - 

Ml-SORT (MX1*MX1*MY1*MY1) 
M2-S0RT(MX2*HX2*rtY2*MY2) 

FM1-M1/F1 
FM2-M2/F 2 

PRINT  50.  LAM0A»M1»F1,FM1#M2,F2,FM2 

FORMAT  ( F lO-r-7.  64^20^1 — - - - 

LAMOA-LAMOA+PI/6. 

K-K*l 

iF(K  .LE.  31  GO  TO  3 
ZETA-ZET A/2 . 

CONTINUE 

END 
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Appendix  B - Force  equation  for  large  amplitude  waviness 


Recalling  that 


f | _L  dh 

f ■ 2 J+  “5  d9  de 

p h 


where 


h = h + z|sin((H-$j)  + z2sin(26+<J2> 


and  p means 


n 02  tt  02  3n  ^2  5h  ^2 

_4'T  e 4"TandT*T  e T'T 


Changing  the  variable  to  h,  another  way  to  calculate  f will  be 


f . liXL. 
r 2 


P h 


. iivjs_  r_Li 

^ “h2  p 


where  p becomes 


h + z'sin(-^r  + \)  -z2<h<h  + z^sin(-£  + \)  + 


and 


h + z|sin(^  + \)  - z2  < h < h + z{sin(-^  + \)  + z2 


'1  4 


where 


\ = $ 

12 


d] 


(2: 


(3 


Consider  only  the  case  zj  <<:  z 2>  where  z2  not  necessarily  small  compared 


with  h,  then  p becomes  approximately 


h-z2<h<h+z2 


and 


h-z2<h<h+z2 


insertion  into  Eq.  (2)  yields 
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where 


when  e « 0,  Eq.  (4)  converges  to 


f 


3 

2uVL  * 

^Z2 


(5) 


which  is  the  case  for  small-perturbation  analysis. 


I 


Appendix  C - Moment  equation  for  large  amplitude  waviness 


Recalling  that 


2hJ  p [ 1+Cesin(0f$j)  + esin(204-02) Y 


where 


For  convenience,  write  M^  as 


^ { INT(M^)  + INT^)) 


where  INT(Mj) , INT(M2)  denote  first  and  second  mode  components  of  the  integral 
in  Eq . (1)  . 

Consider  INT(M^) , when  £ « 1 as  assumed,  it  becomes  approximately 

z ' cos  (0+f  ) sin0 


J»  ^ UT  T.  / O 1 II  V 

-i i r 

p [l+esin(20+$_)  1 


where  p means 


n 02  < o < n 02  h 30  02  < ft  < 5TT  02 

‘4'T  e 4'TandT‘T  e T"T 

by  changing  variable  0 to  9'  as  follows 


6'  - 26  + tf2  + 2 


The  integration  interval  becomes 


0 < 0'  < it  and  2tt  < e'  < 3n 
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■ 


-3b- 


Eq.  (3)  can  be  changed  to  the  following  form 


INT(M1) 


i r 

— j-sin$ 


|,TT  d9  , l'n 

• J»  + eos<*'-'2,i 


cos9d9 

[l-eco89l3 


The  integrals 


+sin($j-02) 
In  Eq . (4)  can  be 


T sln9d6  \ 

I 3 I 

o [l-6cos9l 

carried  out  [ill 


by  writing 


COsf3 


- € + cos9  „ G + cosP 

— or  cos9  » — k- 

1 - G COS9  1+G  COSP 


(4) 


(5) 


such  that 


d0 


(1-e2)  -dp 

(1  + GCOSp) 


results  will  be 


,.TT 

de 

l+(*2/2)  _ 

o (1-Gcos9)^ 

(l-«2)5/2 

,,TT  cos9d6 

(3/2)gtt 

J 3 

O (1-6COS0) 

a-«2>5/2 

|,TT  sin8d9 

2 

J 3 “ 

o (1-GcosO) 

„ 2,2 
(1-e  ) 

hence 


INT(M1) 


, ( TTsln^i  1+(e2/2)  30008(^-02)  e 

’ll'  2 <,.,2,5/2  4 (1.e2)5/2 

sin($1-02)^ 


when  G w 0,  It  converges  to 

INT(Mj)  " *j|-jnin*j  + sin(fj-02)} 


(b) 

(7) 

(8) 


(9) 


(10) 


V 


which  is  the  same  as  derived  for  small-perturbation  analysis. 
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As  for  integral  INT(M2)  , the  C«sin(0+4j)  term  cannot  be  neglected, 
in  fact,  it  is  the  dominant  term  in  the  Integral.  A closed  form  solution 
seems  impossible;  instead,  it  will  be  carried  out  by  series  expansion  as 
follows 

1 


■ 1 - 3y  + 6v2  - 10y3  + 15y^  -+.... 

(i+y)J 

where  in  this  case 

y « e[Csin(&4-$j)  + sin(2&f0^)] 
since  C <<  1,  higher  order  terms  are  approximately 

y2  - e2[2Csin(&t-*1)sin(2e+tf2)  + sln2(29f02)] 

y3  - e3t3Csin(&f*1>sin2(20+02)  + sin3(2fr*?2) ] 
yn  « 6n[nCsin(e+$1)sinn(2&+(?2)  + sin"(20f02) ] 

It  can  be  shown  that,  for  integer  n 

^ + cos(2(M-02)  sinn(2&f02)  sin0d8  « 0 
P 

j + cos(2&f02) sinn(20+02) cos0d9  * 0 
P 

By  changing  variable  as  in  the  evaluation  of  INT(M^),  one  finds 
INT(M-)  - 2z ' ' -3e!Isin(0+\)cos20sin9d9 

L 1 v vl  + 

p 

2 1*  2 

+12e  ; sin(0f\)cos20sin  29sin0d0 
v + 
p 


(11) 


(12) 


-30  c3|  sin(0+^)cos20sin320sin9d8  +- . . .}  (13) 

J *r  > 


where 


X “ *1  * f 


The  final  result,  taking  only  the  first  three  terms,  is 
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INTO^)  “ zj  — 6e[cos*j  - ^cos(*1-02>]  +8  «2sin(*j-02) 
- 20e'£cos*1-Y^cos(*1-<Z>2)  ~\  +-...} 


(14) 


It  is  clear  that  INT(M2>  is  also  a linear  function  of  zj,  and  will  converge 
to  zero  when  € approaches  zero.  From  Eq.  (2,9,14),  the  total  value  will  be 

^ 2h3  ltL  2 3 (1-e2) 3 2 4 (l-.2)5'2 

sin(4  -0  )n  r n 

+ j-j — I + -6e(cos4  ~rcos(4  -0.)) 

(1-e  ) • 

+ 8e2sin($1-02)  - 20e^cos41  -y£c°s(4i-02)) 

Hy  can  be  obtained  in  the  same  way  with  the  following  result 
3„  rr-  ..,  2 


M i JT TT  . 1-K.V2)  3TTsin(42-02) 

“y  " ,r3  zlU_2cos$l  .25/2  + 4 

2h  ( 1 - e ) 


cos($  -0  )-, 

+ 2~2 j + ! -6e(sin*  -^sin(*  -02)) 

(1-e  ) 

+8e2cos(4j-02)  -2Oe3(sin41-|^sin(41-02)) 


When  e « 0,  it  converges  to 


3 

U.VL  R .ITT  , . .1 

^ _3  z^cos$1  + cos(41-02)f 

2h 


(15) 


(16) 


(17) 


which  is  the  case  for  small-perturbation  analysis. 


Appendix  D - Program  for  determining  maximum  thermoelaatic  effecta 


10 


15 

20 


100 


PROGRAM  HEAT (OUTPUT) 

PI»3. 1415926 
DEGRE-57. 296038 
PI2-PI/2. 

0EL--PI2-PI/12-  

BeL-DEUPI/12. 

IFCDEL  .GT.  ( P 12*. 0001 ) ) GO  TO  100 
PHI2»-2.*PI-PI/2A. 

PHI2-PHI2+PT/24. 

TANF1»S IN ( PNI 2) /( COS (PHI2 ) ~-PI2) 
PHIl-ATAN(TANFl) 

SHIFT-PHI2/2.-DEL-PHI1 

PSI-PHI2/2.-DEL 

DELD-OEt*DEGRE 

PHI  20 • PHI  2 *DEGRE 

PHI1D"PHI1*DEGRE 

SHIFO-SHI  FT*PEGRE  _ . . _ _ 

PSID-PSI*DEGRE 

ZWl-PI2*SIN(PHl2/2.  -DEL)  ♦ SI  N ( PH  12  / 2.  ♦ OE  L > 
ZW2  *P I 2-COS (PHI?) 

ZWAVE-ZW1/ZW2  . _ 

FPHAS-C0S(PHI1-PHI2T  ♦ P 12* COS ( PHI 1 ) 
IF(ABSfZWAVE)  .IT.  .0001)  GO  TO  15 
AB-COS(PHl2/2.-OEl)*FPHAS/(ZWAVE*COS (PHIl) ) 
IF( ABS(SIN(PHI1))  .LT.  0.000001)  GO  TO  10 
C-FPHAS/SINCPHIl) 

D-AB*C 

ZPPIM-ZWAVE/FPHAS 

IF ( ZWAVE  .LT.  0.000001)  GO  TO  10 

PRINT  20»DEL0*PHI1D#PSID»ZWAVE»AB*ZPRIM»0  * 

F0RMAT(4F14.A,3E20.4) 

IF ( PHI2  .GT.  (2 , *P I-. 0001 ) ) GO  TO  1 
GO  TO  10 
CONTINUE  "" 

END 


l' 
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Appendix  E - Remarks  on  large  amplitude  wavlness 

If  the  restriction  of  small-perturbation  is  removed,  which  means  first 
mode  and  second  mode  waviness  as  well  as  stator  wave  amplitude  not  necessarily 
small  relative  to  mean  film-thickness,  force  and  moment  will  be  different  but 
not  change  much  if  ^/h  is  not  very  far  away  from  zero. 

Following  the  same  procedure  as  did  before,  the  stator  amplitude  is 


found  to  be 


where 


(-JL)  (I)  (60  2 Ifilfifl.) 

\ 2'  \~M  N ' '■INT(F); 
2tt  w 


P [Ccos(9+$)  + 2cos(20+0  )]cos9d8 

1NT(M)  =*  I T“ 

p [l+Cesin(0+-4j)  + esin^&f^)! 

P [Ccos(0+4)  + 2cos(2&f0„) Id© 

INT(F)  - J -5 

P [l+Cesin(6+i)  + esin^Of^) 1 


We  can  still  conclude  that  Q « 1 in  order  to  keep  within  practical  range 
[Fig.  (9)  and  Appendix  A]. 

Under  the  assumption  of  z^  ^ Zj,  where  may  be  comparable  with  h, 
Appendix  B shows  that  force  will  be  approximately 


2 2 

h3  (1-0 

Therefore  force  will  be  changed  by  the  factor 


/,  2 2 
(1*6  ) 

for  example,  if  c * 0.3,  force  will  be  1.2  times  compared  with  small- 
perturbation  case. 

Appendix  C has  derived  the  moment  integrals.  Results  show  that  moment 
is  a linear  function  of  zj,  and  will  not  change  too  much  if  e is  not  too  big. 
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Appendix  F - Stator  motion  in  Eulerian  Coordinates 

The  formulation  of  Kupperman  (12)  has  been  reviewed,  and  is  felt  to 
be  appropriate  for  small,  steadily  precessing  tilt  of  the  type  studied 
here.  If  one  refers  to  any  advanced  dynamics  text,  one  finds  the  angular 
velocity  vectors,  referred  to  Eulerian  axes,  to  be,  using  present  notation: 


Q = GsinV 
x 


• • 

* Gcosy+0 

Here  6 would  correspond  to  ui  the  speed  of  rotation  of  the  rotor  and  if  the 

stator  follows  rotor  tilt  it  will  be  the  speed  of  precession  of  the  stator. 

The  quantity  0 is  the  spin  of  the  stator  relative  to  the  moving  coordinate 

system  and  must  be  such  as  to  keep  the  absolute  0 = 0,  if  gimbals  or  0-rings 

z 

keep  the  seal  ring  from  rotating.  If  the  tilt  angle  is  fixed  (steady 
precession,  no  rotation),  then  y * 0.  Using  this  information  the  only 
moment  is  that  given  in  Eq.  (10) . 
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